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ABSTRACT 

We have measured the rotational velocities (v sini) of the mass donors in a number of 
D'-type symbiotic stars, using the cross-correlation function method. Four from five 
D' symbiotic stars with known v sini, appeared to be very fast rotators compared 
with the catalogues of v sini for the corresponding spectral types. At least three of 
these stars rotate at a substantial fraction {> 0.5) of the critical velocity. This means 
that at least in D'-type SS the cool components rotate faster than isolated giants. If 
these binary stars are synchronized, their orbital periods should be relatively short 
(4-60 days). We also briefly discuss the possible origin of the rapid rotation and its 
connection with mass loss and dust formation. 
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1 INTRODUCTION 

The Symbiotic stars (SSs - thought to comprise a white 
dwarf (WD) accreting from a cool giant or Mira) represent 
the extremum of the interacting binary star classification. 
They offer a unique and exciting laboratory in which to 
study such important processes as (i) mass loss from cool 
giants and the formation of Planetary Nebulae; (ii) accre- 
tion onto compact objects, (iii) photoionisation and radia- 
tive transfer in gaseous nebulae, and (iv) nonrelativistic jets 
and bipolar outflows (i.e. Kenyon 1986; Corradi et al. 2003). 

Soker (2002) has shown theoretically that the cool com- 
panions in symbiotic systems are likely to rotate much faster 
than isolated cool giants or those in wide binary systems. 
However, there are no systematic investigations of v sini of 
the mass donors in SS. 

On the basis of their IR properties, SS have been clas- 
sified into stellar continuum (S) and dusty (D or D') types. 
The D-type systems contain Mira variables as mass donors. 
The D'-type are characterized by an earlier spectral type 
(F-K) of the cool component, and lower dust temperatures. 
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Among 188 objects in the latest catalogue of symbiotic stars 
(Belczyhski et al. 2000), there are only seven that are classi- 
fied as D'-type: Wrayl5-157, AS 201, V417 Cen, HD 330036, 
AS 269, StHa 190 and Hen 3-1591 (Hen 3-1591 can be S or 
D'). Three of these have been studied using model atmo- 
spheres and all display rapid rotation and s-process elemen- 
tal over- abundances (see Pereira et al. 2005). 

Our aims here are: (1) to measure the projected rota- 
tional velocities {v sini) and the rotational periods (Prot) 
of the giants in D' symbiotic stars, using a cross correlation 
function (CCF) approach; (2) to test the theoretical predic- 
tions that the mass donors in SSs are faster rotators than 
the isolated giants or those in wide binary systems; (3) to 
provide pointers to the determination of binary periods (as- 
suming co-rotation). This is the first of a series of papers 
exploring the rotation velocities of the mass donating (cool) 
components of SSs. 



2 OBSERVATIONS 

The observations were performed with FEROS at the 2.2m 
telescope (ESO, La Silla). FEROS is a fibre-fed Echelle spec- 
trograph, providing a high resolution of A/ A A =48000, a 
wide wavelength coverage from about 4000 A to 8000 A 
in one exposure and a high overall efficiency (Kaufer et al. 
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1999). The 39 orders of the EcheUe spectrum are registered 
on a 2kx4k EEV CCD. AU spectra are reduced using the 
dedicated FEROS data reduction software implemented in 
the ESO-MIDAS system. 



3 « sini MEASUREMENT TECHNIQUE 

Radial velocities and projected rotational velocities have 
been derived by cross-correlating the observed spectra with 
a KO-type numerical mask yielding a cross-correlation func- 
tion (CCF) whose centre gives the radial velocity and whose 
width is related to the broadening mechanisms affecting the 
whole spectra such as stellar rotation and turbulence. De- 
tails of the cross-correlation procedure are given in Melo et 
al. (2001). 

In order to use the width (cr) of the CCF as an estimate 
of V sini one needs to subtract the amount of broadening 
contributing to a unrelated to stellar rotation (e.g., convec- 
tion, instrumental profile, etc.), i.e., (tq. Melo et al. (2001) 
cahbrated cro as a function of [B — V) for FEROS spectra of 
stars with 0.6 < {B — V) < 1.2. The v sinz measured from 
our CCFs for a set of standard stars within this B — V range 
are in good agreement with the literature values. Therefore, 
for all 4 giants in Table^] the Melo et al. (2001) calibration 
has been adopted. 

For V sini greater than about 30 km s^^ the shape of 
the CCF becomes gradually closer to the Gray rotational 
profile (Gray 1976). Therefore in order to correctly fit the 
CCF, a different approach is needed as described in Melo 
(2003). The CCF is then fitted by a family of functions 
CCFv sini = C — D[go®G(V smi)] which is the result of the 
convolution of the CCF of a non-rotating star go, which can 
be fairly approximated by a gaussian, and the Gray (Gray 
1976) rotational profile computed for several rotational ve- 
locities G(V'sini). For each function CCFvsini we found 
the radial velocity Vr, the depth D and the continuum C 
for minimizing the quantity Xvsini which is the traditional 
function with ai = 1, where (Ji is the measurement error 
(see Fig. 1 of Melo 2003, for an example of the procedure). 
The CCFs of our objects are plotted in Fig. 1. 

For V sini> 30 km s~^ the typical error of our v sini 
measurements is 10%. For 10 <v sini^C 30 km s~^ the error 
is about 1.5 km s^^. Our measurements are given in Table 1. 



4 ROTATION OF THE MASS DONORS 

4.1 Individual objects 

HD 330036 : Pereha et al. (2005) obtained L=650 
for the cool component (with possible uncertainties 160< 
L <3000 Lq), Te// =6200±150 K, logg = 2.4±0.7. 
This imply Rg = 22Rq (using L = iTvassRlT^f f), 
Mg =4.46 Mq (using Rg and log g), and Prot < 10.4±2.4 d 
(using Prot V sini< 27ri?g). 

Hen 3-1591 : Medina Tanco & Steiner (1995) give 
spectral type Kl for the cool component. We assume that 
it is luminosity class III. The average radius of a KlIII 
star is 23.9±3 Rq and the average T^// =4280±200 K 
(van Belle et al. 1999). A KlIII star would have a mass 
of 3.9±0.3 Mq (Allen 1973). The uncertainties correspond 



to ±0.5 spectral types. This will resuh in L=172 Lq(±20%) 
and Prot < 51.0±11.3 d. 

StHa 190 : The cool component is of type G4 III/IV 
with T^ff =5300±150 K, logp = 3.0±0.5, and L=45 Lq 
(Smith et al. 2001). This implies Rg = 7.9 ± 0.4 Rq and 
Prot J; 3.8±1.2 days (upper limit calculated for Rg = 
8.3 Rq and sini = 1.0). The upper limit for Prot is con- 
siderably shorter than the supposed orbital periods of 171 d 
(Munari et al. 2001) or 38 d (Smith et al. 2001). 

V417 Cen : Van Winckel et al. (1994) detected a pho- 
tometric period of 245 days. For the cool component they 
obtained G9 Ib-II, log L/Lq =3.5, Teff =5000 K, log 5 = 
1.5±0.5. This implies Rg = 75Rq and Prot < 50.6±10.2 
days. Prot is considerably shorter than the period obtained 
from photometry. However, the photometric period is not 
confirmed with radial velocity measurements and we do not 
know whether this is the orbital period. 

AS 201 : Following Pereira et al. (2005), the cool com- 
ponent is of type F9III with Te// =6000±100 K, L=700 Lq 
(with possible uncertainties 300 < L < 1200 Lq), and 
logg = 2.3±0.3. This imply Rg = 24.5 Rq and Prot < 
49.5±11.0 days. 

4.2 Projected rotational velocities v sinz- 
comparision with catalogues 

There are no systematic investigations of the rotation of the 
mass donors in SSs. The rotational velocities of 13 S-type 
systems listed in Fekel et al. (2003) are between v sinz= 
3.6 - 10.4 km s~^. All D'-type systems so far observed (see 
Table I^J rotate with v sini>20 km s^^. 

The catalog of de Medeiros et al. (2002) of v sin i of lb 
supergiant stars contains 16 objects from spectral type G8- 
KO Ib-II. All they have v sini in the range 1-20 km s"'^. It 
means that V417 Cen is an extreme case of very fast rotation 
for this spectral class. 

The catalogue of rotational velocities for evolved stars 
(de Medeiros et al. 1999) lists ~100 KlIII stars, and 90% of 
them rotate with v sini <8 km s~\ There are only 5 with 
V sin i>20 km s~^ . This means that Hen 3-1591 is a very fast 
rotator (in the top 5%). The same catalog contains 5 objects 
from spectral type F8III-F9III. They rotate with v sini of 
10-35 km s~^. AS 201 is well within in this range. However 
HD 330036 is an extremely fast rotator. The same catalog 
lists 60 objects from spectral type G3,G4,G5 III-IV. They 
all rotate with v sini < 24 km s~^. Again, this means that 
StHa 190 is an extremely fast rotator. 

Thus overall, 4 out of 5 D'-type SSs in our survey are 
very fast rotators. 



4.3 Critical speed of rotation 

There is a natural upper limit for rotation speeds, where 
the centripetal acceleration balances that due to gravita- 
tional at traction, o ften na med th e "critical speed", where 
Vcrit = y/GM/1.5R = 357y/M/R km s'^ (the factor of 1.5 
appears from the assumption that at critical rotation speeds 
the equatorial radius is 1.5 times the polar radius, R). The 
calculated Vait is included in Table |21 No star can rotate 
faster than its critical speed, however we can see that at 
least three D'-type SSs are rotating at a substantial fraction 
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Table 1. Journal of observations and projected rotational velocities of D' type SSs (note that value for AS 201 is from the literature). 
The spectral types of the giants are from different papers (see Section |4. II . (B-V)o is the intrinsic colour for the corresponding spectral 
type. V sini is the rotational velocity measured in this paper. The last column gives other measurements of v sini (if available). 



object 


date-obs 


MJD 


Exp. time 


Cool Star 


(B-V)o 


V sinj 


other 










Spectrum 




[km s-i] 


[km s^-*-] 


HD 330036 


2004-05-22 


53147.212 


2x10 min 


F8III 


0.90 


107.0±10 


lOOitlO" 


Hen 3-1591 


2004-07-19 


53205.222 


2x20 min 


KlIII 


1.09 


23.71±1.5 




StHQ 190 


2004-06-03 


53159.397 


2x10 min 


G4III/IV 


0.88 


105.0±10 


lOOitlO' 


V417 Cen 


2004-04-14 


53109.280 


2x10 min 


G9Ib-II 


0.98 


75.0±7.5 




AS 201 








F9III 


0.58 




25±5" 



"Pereira et al.(2005); ''Smith et al.(2001); 




of their critical speeds: ~ 0.67 (HD 330036), ~ 0.54 

(StHa 190), ~ 0.71 (V417 Cen). Probably for these three SSs 
the orbital inclination is high i 50'' . For the remaining two 
objects we can not exclude the possibility that they also ro- 
tate very fast but are observed at low inclination {i < 30''). 



5 SYNCHRONIZATION AND BINARY 
PERIODS 

5.1 Synchronization in SS 

The physics of tidal synchronization for stars with convec- 
tive envelopes has been analyzed several times (e.g. Zahn 
1977 and also the discussion in Chapter 8 of Tassoul 2000) . 
There are some differences in the analysis of different au- 
thors, leading to varying synchronization time-scales. Here, 
we use the estimate from Zahn (1977, 1989): the synchro- 
nization time-scale in terms of the period is 



: 800 



9 ^ MwD 



+ 1) 



-Port years 



(1) 



where Mg and Mwd are the masses of the giant and white 
dwarf respectively, and Rg and Lg are the radius and lumi- 
nosity of the giant (all in Solar units) . The orbital period P 
is measured in days. 

The S-type SSs are very likely synchronized (Schild et 
al. 2001; Schmutz et al. 1994). Other proof of this suppo- 
sition is that most of the SSs with derived orbital parame- 
ters (see Mikolajewska 2003) have orbital eccentricity e ~0. 
Because the circularization time of the orbits is ~10 times 
longer than the synchronization time (Schmutz et al. 1994 
and the references therein), if a SS's orbit is circularized it 
will very likely be synchronized too. 



A typical D'-type SS would have Rg 20RQ,Lg 
5OOL0 , Alg ~ 3Mq , and white dwarf mass Mwd ~ lAf© . 
For a period of 100 days for a typical D'-type SS, we find 
Tsyn ~ 9 X 10* years. 

For the individual systems, we calculated the synchro- 
nization time (Tsyn) for two derived assuming 
Port =Prot, and Tsyn (100) assumiug Porb =100 days. These 
are given in Table |5| Depending on the individual parame- 
ters, the synchronization time can be from <100 yr up to 
>10'' yr. This means that it is possible for a D'-type SS to 
be synchronized if the orbital period is short (Porb ~Prot). 



5.2 Evolutionary status of the mass donors 

The mass of the mass donors in S-type SSs with known 
parameters are in the range 0.6 - 3.2 Mq (Mikolajewska 
2003). The calculated masses of the mass donors in D'-type 
SSs are larger. As can be seen in Table 2, they range from 
2.2 up to 6.5 Mq. 

Masses of 8 Mq are generally considered the upper limit 
for evolution to planetary nebula nuclei and white dwarfs, 
after heavy mass loss, especially during their AGB phases. 
Following our calculations for the mass of the giants (Ta- 
ble I^J and assuming }Awd ^ 1.4 M©, the total mass of the 
binary is about (Mg-fMi^/u) ~ 3.5 — 8.0Mq, in agreement 
with the above upper limit of 8 Mq for the WD progenitor. 

The position of the mass donors on the H-R diagram 
is presented in Fig|5| assuming near solar chemical com- 
position and stellar parameters as given in Sect 14. II (see 
also Pereira 2005). The evolutionary tracks of Shaller et 
al. (1992) have been used. The donors appeared in a wide 
mass interval - from 2.5 to 7Mq. The derived evolutionary 
masses are in good agreement with those obtained from Rg 
and fog 3. Three of them appeared crossing the Hertzsprung 
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Table 2. Parameters of D'-type SS. IR types are from the catalogue of Belczynski et al.(2000), v sin i is the rotational velocity of the mass 
donors (as adopted here), Rg and Mg are the radius and the mass of the giant, respectively (see Sect 14. ll . Tsyn is the synchronization 
time of the binary adopting Porb =Prot, Tsyn(lOO) is the synchronization time of the binary adopting Po^b =100 days, Vcrit is the critical 
rotational velocity of the giant, Prot is the rotational period of the giant, a is the semimajor axis of the system calculated supposing 
synchronized rotation {Porb = Prot)- The values in brackets for Prot and a are the upper limits, supposing 10% uncertainity in v sini 
and Rg. 



object IR Cool Star V sini JJ^rit R-s Mg Tsyn Tsj/n (100) Prot a, 

type Spectrum [km s"-*-] km s~^ [R©] Mq [yr] [yr] [days] [R-©] 



HD 330036 


D' 


F8III 


107.0 


160 


Hen 3-1591 


S,D' 


Kim 


23.7 


144 


StHa 190 


D' 


G4III/IV 


105.0 


191 


V417 Cen 


D' 


G9Ib-II 


75.0 


105 


AS 201 


D' 


F9III 


25.0 


150 



22.1 


4.46 


25 


2.1. 10^ 


<10.4 (12.8) 


35(40) 


23.9 


3.9 


8680 


1.3. 10^ 


<51.0 (62.3) 


98(112) 


7.88 


2.25 


27 


1.3.10'^ 


<3.78 (4.64) 


15(17) 


75.0 


6.45 


36 


553 


<50.6 (61.8) 


112(128) 


24.5 


4.35 


6423 


1.1. lO'"^ 


<49.5 (60.5) 


99(113) 



5.3 Clues regarding the orbital periods 



-6.0 




3.4C' 



Figure 2. The position of the mass donors on the H-R diagram 
(see Sect l5.2t . From top to bottom, the objects are V417 Cen, 
AS 201, HD 330036, Hen 3-1591, StHa 190. The typical error 
on the bolometric magnitude, M-bol^ is plotted in the bottom left 
corner. The typical error on logT^^j is ±0.01. 



gap (HD 330036, V417 Cen and AS 201), StHa 190 is situ- 
ated near the base of the red giant branch, and He 3-1591 
is already evolving on the red giant branch. 

The relevant time for a 5M0 star to cross the 
Hertzsprung gap is ~8xl0^ yr and its life time on the red 
giant branch is RiSxlO^ yr (Iben 1991). For a 1.5Mq star 
these times are 1.5x10* yr and 1.57x10* yr (Iben 1991). 
These lifetimes are longer than the calculated Tsyn and com- 
parable with Tsyn (100). This means that the rotation of the 
mass donors in D'-type SSs could be synchronized for these 
lifetimes. 



Because the orbital periods of the majority of SSs are 
unknown, an indirect method to obtain Porb is to mea- 
sure V sini. If the mass donors in SSs are co-rotating 
(Prot ^Porb), wc cau fiud clues for the orbital periods via 
the simple relation Porb Vrot = 27rRg, where Port is the or- 
bital period, Vrot and Rg are the rotational velocity and 
radius of the giant, respectively. It is very useful in the case 
of eclipsing binaries, where sini « 1 and v sinifsvrot- If the 
inclination is unknown, we can only obtain an upper limit 

for Vorb- 

Using Kepler's third law [47ra^ = G{Mg + Mwd)P^], 
we can calculate the semimajor axis of the systems. These 
are given in Table |21 The values in the brackets (for Prot 
and a) correspond to the estimation of the upper limit of 
these parameters, assuming 10% errors in v sini and Rg. 

Up to now, from 188 SSs, the orbital elements and bi- 
nary periods are well known for ~40 objects only (and they 
are all S-type). The orbital periods are in the range 200 - 
2000 days (Mikolajewska 2003). As can be seen in Tabled 
if D'-type SSs are synchronized, their orbital periods would 
be relatively short (4-60 days) and the distance between the 
WD and the mass donor would be 2-5 Rg. 



6 DISCUSSION 

The observation of fast rotation of D'-type SSs raises two 
further questions; first, what is the evolutionary history of 
these stars which has produced such high rotation, and sec- 
ond, what does the effect of high rotation have on their mass 
loss and subsequent dust formation? 

We see three possible reasons for the fast rotation of 
mass donors in D'-type SSs: 

(i) the rotation is synchronized with the binary period 
(Prot = Porb)- In this case their orbital periods should be 
short <50 days. 

However, it also could be, that they are not syn- 
chronously rotating. The possibility that Prot > Porb has 
to be excluded because the orbital separations would be un- 
reasonably small, and the synchronization time would be 
extremely short (they will be synchronized in 30-9000 years; 
see Tsyn in Table |5J. If Prot < Porb, reasons for their fast 
rotation could be: 
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(ii) the current giants have been spun-up from the 
transfer of angular momentum. Jeffries & Stevens (1996) 
proposed a mechanism in which accretion of a slow mas- 
sive wind from the AGB progenitor of the current white 
dwarf can transfer sufficient angular momentum. This also 
explains the chemical enrichment in s-process elements in 
D' SSs, that were present in the AGB wind (see also Pereira 
et al. 2005). Mass transfer via Li, when the current white 
dwarf was the mass donor, can also spin-up the companion, 
as in millisecond radio pulsars (van den Heuvel 1984). 

(iii) planet swallowing to spin up the giant. Rough esti- 
mation gives angular momentum transfer during a collision 
of a planet with mass nip and velocity Vp to a giant of mass 
Mg, of the order of AF = rUpVpRg. This causes a change 
in the giant's rotational velocity Au^"' — rUpVp / const. Mg, 
where cansX.=Jg/MgR^ depends on the internal structure of 
the giant, [const. = 0.4 for a solid uniform density sphere, 
and less for a star-like centrally condensed sphere). Assum- 
ing a centrally condensed star (giant) such that const. = 
0.01, Vp = ^GMg/R ~ 10 - 100 km s-\ A/g = 2 - lOA/® 
and nip = O.OIM© we estimate Au^"' ~ 1 — 50 km s"^, 
showing that in the right circumstances, the planet could 
spin up the giant to the rotational velocities observed in 
D'-type SSs. 

Fast rotation, i.e. Vrot ^ 0.5wcrit, may change a spheri- 
cal star with a spherical wind into an equatorially flattened 
system, with both the radius of the star and stellar wind pa- 
rameters depending on the stellar latitude. Such stars will 
have an equatorial radius significantly larger than the polar 
one, and equatorially enhanced mass loss (see Lamers 2004 
and references therein). 

Because it seems that the majority of the giants in 
D'-type SSs are rapid rotators fSect. H^ . we expect that: 

(1) they have a larger mass loss rate than the slower 
rotating giants; 

(2) their mass loss is enhanced in the equatorial regions; 

(3) they could be even equatorially flattened. 

It is possible, that the dusty environment in D'-type SSs 
is connected with rapid rotation of the mass donors. Intense 
mass loss in the equatorial regions can lead to the formation 
of an excretion disc in which the higher gas density enhances 
dust formation and growth. The broad IR excess in D'-type 
SSs can be due to the temperature stratification in the dust 
from such an excretion disc (Van Winckel et al. 1994). Other 
possible explanations for the presence of dust can be that 
it is left over from the formation of a planetary system, or 
it is a relic from a strong dusty mass loss when the present 
day white dwarf was on the AGB. However, we consider it is 
more likely that it originates in the current outflow and that 
this is enhanced in the equatorial regions by rapid rotation. 



7 CONCLUSIONS 

Our main results are: 

(1) We have measured the rotational velocities of the 
mass donors in D'-type symbiotic stars, using the CCF ap- 
proach. 

(2) Four of the five objects appeared to be very fast 
rotators compared with the catalogues of v sin i for the cor- 



responding spectral types. At least three of them rotate at 
a substantial fraction (> 0.5) of the critical velocity. This 
means that at least in D'-type SSs, the cool components ro- 
tate faster than isolated giants (as predicted by Soker 2002). 

(3) If they are tidally locked, the orbital periods should 
be as short as < 10-50 days. 

(4) As a result of the rapid rotation, they must have 
larger mass loss rates than the more slowly rotating giants, 
and their mass loss is probably enhanced in the equatorial 
regions. 

To understand better these objects, we need their bi- 
nary periods to be derived from radial velocity measure- 
ments and the inclination determined. In subsequent pa- 
pers, we plan to explore the projected rotational velocity of 
the cool giants in the other types of symbiotic stars and to 
compare their rotational velocities with that of the isolated 
giants with similar mass and evolutionary stage. 
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